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Abstract 

Control by dissipation, or environment engineering, constitutes an important methodology within quantum coherent control 
which was proposed to improve the robustness and scalability of quantum control systems. The system-environment coupling, 
often considered to be detrimental to quantum coherence, also provides the means to steer the system to desired states. This 
paper aims to develop the theory for engineering of the dissipation, based on a ground-state Lyapunov stability analysis of 
open quantum systems via a Heisenberg-picture approach. In particular, Lyapunov stability conditions expressed as operator 
inequalities allow a purely algebraic treatment of the environment engineering problem, which facilitates the integration of 
quantum components into a large-scale quantum system and draws an explicit connection to the classical theory of vector 
Lyapunov functions and decomposition-aggregation methods for control of complex systems. This leads to tractable algebraic 
conditions concerning the ground-state stability and scalability of quantum systems. The implications of the results in relation 
to dissipative quantum computing and state engineering are also discussed in this paper. 
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1 Introduction 

Control of quantum systems lies at the core of the 
quantum technology [46,6,1], while stability analysis 
provides the appropriate tool for the systematic de¬ 
velopment of quantum control theory. The stability 
analysis has been used in several quantum control syn¬ 
thesis problems [24,9,31,22,32,47]. The applications 
include measurement-based feedback control and co¬ 
herent control for the generation of quantum states as 
well as the regulation of system performance. Among 
all the methods for stability analysis, the Lyapunov 
stability approach is the most fundamental, as the en¬ 
ergy of a quantum system is well-defined for most of 
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the physical systems and a Lyapunov function can be 
easily constructed [18,45,35,41,2]. In particular, as we 
will demonstrate in this paper, the Lyapunov method 
provides a means for engineering the dissipation to be 
used as coherent control. 

Quantum computing often involves the execution of a se¬ 
quence of unitary operations on quantum systems. How¬ 
ever, the severe decoherence associated with the quan¬ 
tum systems presents a major obstacle to the scalabil¬ 
ity of this approach. For this reason, methods for ro¬ 
bust realization of unitary operations are currently un¬ 
der discussion. The possible plans include topological 
quantum computing, adiabatic quantum computing and 
dissipative quantum computing. Among these schemes, 
the adiabatic quantum computing and dissipative quan¬ 
tum computing have direct relevance to the stability of 
quantum systems. For example, in dissipative quantum 
computing and state engineering, dissipation is intro¬ 
duced as a resource to coherently control the system [44]. 
The idea is to consider open quantum systems, and sta¬ 
bilize their quantum states by engineering the system- 
environment interaction. If designed judiciously, the dis- 
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sipation will drive the system to a target steady state re¬ 
gardless of the initial state and external perturbations. 
This method can be used to generate highly entangled 
quantum states, and perform quantum computation by 
encoding the computation result to the steady state of 
the system. Since dissipation of energy is the key phys¬ 
ical principle behind this method, this kind of coher¬ 
ent control approach can be referred to as control by 
dissipation. Our goal in this paper is to formulate the 
method of control by dissipation within the framework 
of ground-state stability, and then propose approaches 
for the synthesis of system-environment dissipative in¬ 
teractions that rely on Lyapunov methods for stability 
analysis. 

Stability of quantum states has been the focus of many 
theoretical studies. Many of them have successfully de¬ 
rived sufficient conditions for convergence of quantum 
Markov systems to a steady state [40,7,16,37,15,34,27]. 
In particular, the stability of quantum states in a dis¬ 
sipative setting has been considered in [42,37,43,1]. In 
these studies, the target state is often explicitly given 
and follows a Schrodinger-picture master equation. The 
dissipative couplings, compensated by Hamiltonian con¬ 
trol, can generate a Markov process that converges to the 
target states [42,43]. The implementation of the system- 
environment couplings with the practical resources has 
been investigated experimentally. Dissipative engineer¬ 
ing of several types of quantum systems has been demon¬ 
strated in recent years [3,17,11,20,36,38]. 

In this paper, we adopt an alternative path to approach 
the stability theory within the Heisenberg picture, where 
instead of designating target states explicitly, they are 
characterized as ground states a Lyapunov operator, and 
the stability problem is transformed to the problem of 
stabilization of the ground states of the Lyapunov oper¬ 
ator. The formalism of Lyapunov stability can thus be 
conveniently introduced to engineer the desired system 
dissipation within this framework. This allows to derive 
tractable sufficient conditions expressed in terms of op¬ 
erator inequalities, which can be used for the synthesis 
of the desired system-environment coupling. Such con¬ 
ditions is the main contribution of this paper compared 
to our previous work [27]. The general results regarding 
stability of Lyapunov operators obtained in [27] do not 
readily apply to the problem of control by dissipation. 

An important advantage of the Heisenberg-picture ap¬ 
proach developed here is that the target state does 
not need to be given in advance. In addition to the 
entangled-state engineering applications in which the 
Lyapunov operator is chosen based on the knowledge of 
the target state, there exists a large class of applications 
where the control goals are posed as minimization of the 
expectation of an operator while the target state with 
respect to which the expectation is taken is not known. 
For example, the problems of sequential quantum com¬ 
putation and the quantum satisfiability problem (SAT) 


[5,25] involve operators which play the role of cost 
functions. In these problems, the target states which 
minimize the expectation of the operators are unknown 
and result from computation and/or control. More¬ 
over, the target state in these applications may be not 
unique. This complicates the analysis based on the con¬ 
ventional Schrodinger-picture approach. Therefore, the 
Heisenberg-picture approach extends the applicability 
of the control by dissipation. 

One of the main contributions of this paper is con¬ 
cerned with the scalability of the control by dissipation, 
when this control method is applied to large quantum 
systems comprised of multiple interacting subsystems 
coupled with the environment. The Heisenberg-picture 
Lyapunov approach has an advantage in that the 
problem can be treated in a way that resembles the 
decomposition-aggregation engineering [4,39] for com¬ 
plex classical systems. Namely, a large-scale quantum 
system is decomposed into subsystems and an individ¬ 
ual Lyapunov operator is associated with a subsystem. 
This allows us to establish conditions, expressed in 
terms of the subsystems’ Lyapunov operators, under 
which the quantum system is guaranteed to converge 
to its ground state. Here we note a similarity with the 
classical connective stability conditions [39], which have 
proved to be useful in the synthesis of decentralized 
controllers for large-scale systems. 

A typical methodology for the synthesis of dissipations 
involves two problems, the calculation of the stabilizing 
system-environment couplings and the implementation 
of these couplings using the available physical resources. 
For example, it is possible to construct a coherent opti¬ 
cal network to realize a linear coupling [26]. Therefore, in 
this paper we focus on the first problem of calculation of 
coupling operators that render the states of the ground 
energy asymptotically stable. Particularly, we can apply 
this method to check the feasibility of the solutions pro¬ 
posed in [44]. It is worth mentioning that the constraints 
on the system-environment couplings could be greatly 
relaxed if Hamiltonian control is available [42,43]. 

The preliminary version of this paper has been ac¬ 
cepted for presentation at the American Control Con¬ 
ference [28]. Compared to the preliminary conference 
version, this paper has been substantially revised and 
expanded. It includes a detailed exposition of the back¬ 
ground on open quantum systems, the new material 
on the scalability of the Lyapunov methods, synthesis 
of dissipation, the examples and an application to sta¬ 
bilization of quantum states associated with quantum 
toric codes. The paper also includes detailed proofs of 
all the results and gives detailed discussions of these 
results, which were not included in [28]. 

The paper is organised as follows. In Section 2, we in¬ 
troduce the notations and the model considered in this 
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paper. In Section 3 we present the ground-state stabil¬ 
ity analysis of Lyapunov operators. Section 4 discusses 
the scalability problem, where a large quantum system 
may be governed by more than one Lyapunov operators. 
Section 5 concerns with the synthesis of the dissipation. 
More explicitly, this section concerns with the calcula¬ 
tion of the correct coherent couplings for the ground- 
state stabilization when the Lyapunov operator is given. 
Conclusion is given in Section 6. The proofs of the re¬ 
sults are given in the Appendix in Section 7. 


2 Notations and preliminaries 

2.1 Open quantum systems 

Consider a Hilbert space B and define B{B) as the space 
of bounded operators on B. We only consider finite- 
dimensional quantum systems throughout this paper. 
In other words, B is assumed to be finite-dimensional 
throughout the paper. Hence all bounded operators in 
our case are representable as complex matrices. Let X € 
B(B). X T denotes the transpose of X and is the 
adjoint of X. An operator X is called an observable if 
= X. The notation X > 0 (X < 0) means the opera¬ 
tor A is a positive (negative) semidefinite operator. We 
write X > 0 if X is positive definite. Also, we will use 
the notation X y 0 for positive semidefinite operators 
X whose smallest eigenvalue is equal to 0. 

Given a bounded observable X € B(B ) and a trace class 
operator p on B , ( X) p denotes the trace of Xp, ( X) p = 
Tr Xp. When p is a density state, i.e,, a matrix whose 
trace is equal to 1, then ( X) p is the mean value of X 
evaluated at the density state p. 

Control by dissipation is implemented by coupling the 
systems to a collection of environments. To describe 
evolution of the quantum system subject to an envi¬ 
ronment, in addition to the Hilbert space B, consider 
the environment on a Fock space Bb over L 2 (M + ,<if) 
corresponding to Boson field modes. The system evolu¬ 
tion in the Heisenberg picture is captured through the 
time evolution of observables (self-adjoint operators) of 
the system. More precisely, the observable X evolves as 
X{t) = U(t)^{X®I)U(t), where U{t) is the unitary evo¬ 
lution operator of the combined system. The dynamical 
equation for X(t) can be expressed using the quantum 
stochastic differential equation [29] 

K 

dX{t) = {-i[X, H] + Z(X))dt + X(t)]dB k (t) 

1 

K 

+ ^[X(t), L k (t)}dBl(t), (1) 

fc=l 


with 

K 1 1 
£(X) = £ L\XL k - -L{L k X - -XL\L k . (2) 

fc=l 

Here H is the Hamiltonian of the system, L k describes 
the coupling between the system and the k -th of the total 
K environment fields, B k {-) and b\{-) are the annihila¬ 
tion and creation processes defined on Bb- Equation (1) 
defines a Markov process. The generator of this Markov 
process is determined by 

g(X) = -i[X,H]+£(X). (3) 

In conjunction with the Heisenberg picture dynamics, 
the evolution of the density state p t in the Schrodinger 
picture is given by 

pt = —i[H,p t ]+^^L k p t L\. — -L\L k p t — -p t L\L k . (4) 

k 


2.2 System decomposition 

In the sequel, we will make use of the decomposition of 
the Hilbert spaces B into the tensor product of Hilbert 
spaces B = (ff), B t . Each {Bi} can be thought of as 
a Hilbert space on which subsystem i is defined. As a 
simple example, a system consisting of M interacting 
two-level systems (qubits), can be defined in terms of the 

Hilbert space B = (8>)=i Bi, where Bi = C 2 , the two- 
dimensional Euclidean complex space. This motivates 
considering the mentioned decomposition of B. 

Furthermore, the system observables W \, A = l,2,...,iV, 
which could be associated to a subset of {Bi}, can be 
defined through a standard embedding. For example, if 
some W\ is defined on I (g)... (g) Bi (g>... (g> Bj (8> • • • <S> /, 
then W\ is said to be associated to Bi and Bj, or as¬ 
sociated to subsystems i and j. The coupling operators 
{L k }, k = 1,2, ...,K are defined on B as well. We say 
L k is associated to W\ if [ L k , W\] ^ 0. 

2.3 Lyapunov operators and ground state stability of 
complex quantum systems 

We recall the definition of the Lyapunov operator [27]: 

Definition 1 A quantum Lyapunov operator V is an 
observable (a self-adjoint operator) on a Hilbei't space B 
for which the following properties hold: 

(i) VhO- 

(ii) G{V) < 0. 
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One natural choice of the Lyapunov operator is the en¬ 
ergy operator of the system. For example, the Lyapunov 
operator can be defined by offsetting a system Hamil¬ 
tonian H as V = H — d >: 0, where d is the smallest 
eigenvalue of hB 

In this paper, we restrict our attention to considering 
the observables which satisfy the condition 

[X,H] = 0. (5) 

This condition holds, for example, in the state engineer¬ 
ing by dissipation problems concerned with preparation 
of the ground states of Hamiltonians, when no additional 
Hamiltonian control is used [43,44,30]. In these prob¬ 
lems the observable of interest is X = H — d and (5) 
holds trivially. More generally, under condition (5) the 
evolution of the observable X described by (1) is due to 
the environment. This allows us to focus entirely on the 
analysis and synthesis of the effects associated with the 
environment, which is the main objective of this paper. 

For any observable X which satisfies (5), the expression 
(3) for the system generator is simplified into 

g(X) = ]T L\XL k - l -L\L k X - \xL\L k . (6) 

k 

We use po to denote the initial density state of the sys¬ 
tem, and pt to denote the system state at time t. Fol¬ 
lowing Meyer [23], any convergence of a trajectory in 
the form of pt —> Poo should be understood as conver¬ 
gence of probability distributions, i.e., pt —» /?oo means 
TV p t X —> TV Poo X for all bounded X. This is a nice con¬ 
vergence property which can be studied using the tools 
of quantum probability such as the Prokhorov theorem 
and tightness [23]. Recall that a sequence (p„) n >i of 
trace-class operators on a Hilbert space is tight if for any 
e > 0, there exists a finite rank projection P and no > 0 
such that TV p n P > 1 — e for all n > n$. In addition, 
the Prokhorov theorem states that (p n )n>i is sequen¬ 
tially compact if (p„) n >i is tight. In this paper, since 
H. is finite-dimensional, all projections are finite rank. 
Therefore, any trajectory p t is tight [23,27], and hence 
admits a subsequence converging to a limit point p'. 

The ground-state stability of an operator X is defined 
using the mean of the operator: 

Definition 2 Suppose the smallest eigenvalue of an ob¬ 
servable X is d. X is said to be asymptotically ground- 
state stable if 

(A (t))p 0 = {X) Pt —> d, as t —> oo. (7) 


1 In accordance with the common convention of quantum 
physics, the identity operator is omitted here and elsewhere, 
i.e., H — d should be understood as H — dl. 


Here {X) Pt is an alternative representation of (X(t)) Po 
in terms of pt . 

Consequently, a Lyapunov operator V is asymptotically 
ground-state stable if 

(V(t )) p o = (V) pt -»> 0. (8) 

Remark 3 We note a parallel between the ground-state 
stability and convergence of probability distributions of 
classical Markov stochastic processes to stationary dis¬ 
tributions. Precisely, let x{t) be an R l valued Markov 
stochastic process with the transition probability func¬ 
tion P(t,x,Y) = Prob(x(t) £ T|a;(0) = x), here Y is 
a measurable set. It is known that under certain condi¬ 
tions [14, Theorem 4.3 and p. 121], there exists a unique 
stationary probability distribution p(-) for x(t) such that 
for any bounded function /( x), the expectation of f(x) 
with respect to P(t , x, •) converges to that with respect 
to p: 

lim / f(y)P(t, x, dy) = f f(y)p(dy). 

‘-*- 00 Jr ‘ Jr 1 

This property reads that p is an attracting equilibrium 
distribution for {P(f, x, •), t > 0}. It is analogous to the 
property (X) Pt —>• d = ( X) p > considered above. 

If, in addition, Y 0 is a set with the property P(t, x, I'd) —► 
0, then letting f{x) be the indicator function of Y 0 , 
f(x) = Xy 0 (x), leads to lim^^ f Rl XY 0 (y)P{t,x,dy) = 
f R i XY 0 {y)p(dy) = 0. This property is analogous to 
(8); it further illustrates the analogy between station¬ 
ary equilibrium distributions of Markov processes and 
the ground states of asymptotically ground-state stable 
Lyapunov observables. 

Definition 4 The state trajectory p t is said to converge 
to a set S if the limit points of pt are all contained in S. 

This definition is often used to characterize the con¬ 
vergence to an invariant set in the invariance principle 
[24,45,27]. Denote Zx = {p ■ ( X) p = d} to be the set of 
the ground states of X. 

Proposition 5 The state trajectory p t is converging to 
Zx if and only if X is asymptotically ground-state stable. 

The proof of this and subsequent results are given in the 
Appendix in Section 7. 

We will also exploit the notion of dissipation functional: 

Definition 6 ([21,7]) The dissipation functional of an 
operator X is defined as 

£>(X) = Q(X f A) - g(X^)X - X^g(X). (9) 
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The dissipation functional characterizes the dissipation 
of energy. With a single coupling operator L, the dissi¬ 
pation functional is calculated to be 

Q(X) = [L',X*][X,L], (10) 

and hence 2) (A) > 0. 

The first objective of this paper is to develop the ground- 
state stability theory using the notion of Lyapunov op¬ 
erator. We will consider Lyapunov operators which are 
associated either with a single quantum system, or a sub¬ 
system of the total system. The need in such theory can 
be illustrated by the following result from [27]. 

Proposition 7 [27, Theorem 8] Suppose V is a Lya¬ 
punov operator of the system. The state trajectory p t 
will converge to Zy if ('S(V)) P > 0 for p Zy and 
[G(V),V] = 0. 

As one may see, the conditions in Proposition 7 are 
not fully algebraic and also they are not easy to verify; 
the validation of these condition requires computing the 
mean (T>(V)) P , which in tern requires p to be explicitly 
given or known. As discussed, this may not be possi¬ 
ble in certain applications. This motivates us to revisit 
the ground-state stability theory for Lyapunov opera¬ 
tors, in order to derive ground-state stability conditions 
expressed purely in terms of operator inequalities. 

The second objective of this paper is to apply this the¬ 
ory to stabilization of large-scale quantum systems. We 
consider an operator W = XaLi W'A■ W\ y 0, as a can¬ 
didate for the Lyapunov operator of a large-scale system, 
where each W\ is an observable associated with one of 
the subsystems or a group of the subsystems of the sys¬ 
tem. Such operator sum representations naturally arise 
in many problems of control by dissipation, including 
the preparation of multipartite entangled states, hence 
the theory developed in this paper is aimed at these ap¬ 
plications. For more information about the applications 
of control by dissipation please refer to [44,30,43]. 

In general, the fact that the individual observables W\ 
have zero eigenvalue (as implied by the notation W\ y 0 

VA) does not guarantee that W = XaLi Wa has a zero 
eigenvalue; in fact W can be positive definite. In the light 
of the definition of the Lyapunov operator, this means 
that the operator sum W = XaLi W''a may result in an 
operator which does not satisfy formally all the proper¬ 
ties of Definition 1 (recall that by definition, Lyapunov 
operators have a zero eigenvalue). For this reason, we use 
the notation W, W\ instead of V in Section 4 because 
there is a possibility that W, W\ cannot be made Lya¬ 
punov operators even if they are asymptotically ground- 
state stable. Similarly, for W to satisfy the condition 
Q{W) < 0, the condition Q(W\) < 0 does not have to be 
satisfied for all W\ ; that is, W\ may not be a Lyapunov 


operator either (for W\ to be a Lyapunov operator, the 
condition G(W\) < 0 must be satisfied). We will show 
in Section 4 that the approaches to the engineering of 
the ground-state stability of W can be quite different 
depending on whether or not II'a can be taken to be a 
Lyapunov operator of the subsystems. 

It is worth mentioning that the issues discussed above 
are similar to those arising within the vector Lyapunov 
function approach [4]. In certain situations arising in 
the classical stability theory for large-scale systems, it 
is more convenient to use a vector Lyapunov function 
rather than a scalar function for a large-scale system 
[39]. Indeed, in general, scalar functions comprising the 
vector Lyapunov function of a stable large-scale system 
do not need to be Lyapunov functions individually. Par¬ 
ticularly, the decomposition-aggregation method used in 
[39] to simplify the analysis by decomposing the large 
system into several subsystems made extensive use of the 
vector Lyapunov function machinery. When the subsys¬ 
tems are coupled together, a connective stability condi¬ 
tion will ensure the total system is stable after the ag¬ 
gregation. In our case, W is the quantum counterpart 
of the vector Lyapunov function, and W\ is the quan¬ 
tum counterpart of the scalar component of that func¬ 
tion. Also, the scalability property discussed in this pa¬ 
per where the ground-state stability of the operator IF is 
derived from the ground-state stability properties of the 
addends {IFa}, is parallel to the classical decomposition- 
aggregation approach mentioned above. Since each Wa 
may act on several subsystems and some coupling may 
affect more than one subsystems, a scalability condition 
is needed to ensure the cross-couplings do not undermine 
the stability of W. 


The following lemma summarizes the approaches to the 
stability of a large-scale quantum system. 

Lemma 8 Given a collection of observables IF \ y 
0, A = 1,2, ...,1V, consider W = XaLi^A; whose 
smallest eigenvalue is d. 

(i) Suppose (W\) Pt —> 0 for each A. Then (W) Pt —► 0 
and d = 0. 

(ii) Conversely, suppose (W) Pt —► 0. Then d = 0 and 
each W\ is asymptotically ground-state stable. 

The proof of the lemma is deferred to the Appendix. 


Remark 9 If (W) Pt —> d > 0, each W\ is not neces¬ 
sarily asymptotically ground-state stable. As a simple 
example to illustrate this, define W = W\ + W 2 with 


Wi 



y o, w 2 



y 0. The smallest 


eigenvalue of IF is 1. Letting p t 




we ob¬ 


tain that the condition (W) Pt —> 1 is equivalent to p ] 1 + 
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p 22 —> 1. Now suppose both W\ and W 2 converge to their 
ground states, then we must have (Wi) Pt = p} 1 —> 0 and 
(W 2 ) Pt = p 22 — > 0 simultaneously, thus p) 1 + p 22 —> 0, 
which contradicts the condition (W) Pt —► 1; also see 
statement (i) of Lemma 8. The established contradiction 
shows that W± and W 2 cannot converge to their ground 
states simultaneously, and at least one of them is not 
ground-state stable. 

Lemma 8 suggests two different approaches to engineer¬ 
ing of the ground-state stability of W (as we mentioned, 
investigation of such approaches is the main objective 
of this paper), namely, through engineering the ground- 
state stability of every W\, or the ground-state stabil¬ 
ity of X] AeA ' W\. Here A is a subset of the set of all 
A, A. The trivial case where A is divided into {A, 0} 
means that we engineer the ground-state stability of W 
directly. In this paper by engineering we mean the syn¬ 
thesis of coupling operators between the environment 
and the system. 

We would like to mention that Lemma 8 has a con¬ 
nection with the notion of frustration-free Hamiltonian 
[44]. A Hamiltonian H in the form of H = Y1\H\ 
is called frustration-free if the ground states of H are 
also the ground states of every H\. Suppose H\ y 0, 
d is the smallest eigenvalue of H and p g is one of the 
ground states of H. Then if d = 0, we have (H) p = 
Y^\(H\)p g = 0 and so (H\) Pg = 0, which proves that 
H is frustration-free. Therefore, if W and W\ denote 
Hamiltonians in Lemma 8, the condition (W\) Pi —> 0 of 
(i) and the condition (W) Pt —> 0 of (ii) in fact imply the 
frustration-freeness of W through proving d = 0. 

The property of the system observables to maintain their 
smallest eigenvalue to be d = 0 while adding the sub¬ 
system Hamiltonians H\ and associated observables W\ 
means that the system size can be increased without per¬ 
turbing the ground energy. Such a scalability property 
is often desired in quantum engineering. 

3 Lyapunov stability of the ground states 

In this section, we consider the generator with one dis¬ 
sipation channel 

g{x) = tfxL - ^lx - \xi)l. ( 11 ) 


Recall that a state pi is an invariant state of the quantum 
system, if it satisfies the condition (X(t)) Pl = (X) PI for 
any operator X [7,27]. Thus we have {G{X)) Pl = 0 for 
an invariant state pi. 

The next statement gives the quantum version of the 
Lyapunov’s second method for stability. 


Lemma 10 Suppose V is a Lyapunov operator of the 
system. If {Q{V)) P < 0 for any p Zy, thenV is asymp¬ 
totically ground-state stable. 

The proof of the lemma is deferred to the Appendix. A 
special case of Lemma 10 is concerned with the generator 
satisfying the condition 

Q{V) < -cV, c>0. (12) 

In this case, we can integrate (12) to obtain (V(t)) Po < 
e~ ct (V) Po . The system exponentially converges to the 
ground states of V. 

The exponential convergence condition (12) does not de¬ 
scribe all the dynamics that lead to the asymptotical 
stability of the ground states. Not all physical systems 
are exponentially stable. A more general treatment will 
involve dealing with the condition Q(V) < 0. To this end 
we will make use of the dissipation functional S?(V). 

Lemma 11 IfV is a Lyapunov operator of the system 
satisfying 55 (V) > cV 2 for some c > 0, then V is asymp¬ 
totically ground-state stable. 

Particularly, we can make use of Lemma 11 to obtain 
the following result. 

Lemma 12 IfV is a Lyapunov operator of the system 
satisfying cV < S(V) for some c > 0, then the state 
trajectory pt will converge to Zy. 

With multiple dissipation channels, the generator of the 
Lyapunov operator V is expressed as (6) and the dissi¬ 
pation functional becomes 

nV) = J2l L lVW,L k \. (13) 

k 

All the above stability results can be routinely extended 
to the multi-channel case. 

Now we define two conditions on an observable X as: 

Condition ES An observable X is said to satisfy con¬ 
dition ES, if there exists a constant c > 0 such that 

G(X) < -cX, c > 0. (14) 


Condition DS An observable X is said to satisfy con¬ 
dition DS, if there exists a constant c > 0 such that 

G(X) < 0, (15a) 

cX<D(X),c> 0. (15b) 
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Using these conditions, the sufficient conditions for the 
convergence to the ground states of general quantum 
systems can be expressed by operator inequalities, as 
summarized in the following theorem: 

Theorem 13 IfV is a Lyapunov operator of the system 
satisfying either Condition ES or Condition DS, then the 
state trajectory pt will converge to Zy. 

It is worth mentioning that under Condition ES, The¬ 
orem 13 guarantees that the state trajectory converges 
exponentially, hence the name ES - exponential stabil¬ 
ity. Likewise, under Condition DS Theorem 13 guaran¬ 
tees that the convergence is asymptotic and is ensured 
by dissipation properties of the system, hence the nota¬ 
tion DS- ‘dissipative stability’. 

Below are two examples to illustrate the use of Condi¬ 
tions ES and DS. 

Example 14 Consider a two-level quantum system 


with the system Hamiltonian H = 


. In this 


example, we wish to engineer a dissipative coupling to 
govern the system to the ground state of the operator 
V which is defined as 


Clearly, V = H +1/2 and we have [ V , H] = 0 as required. 

To verify condition DS, consider a general coupling op- 

( Zoo ^oi \ 

and sub- 

Uo hi ) 

stitute this expression in (15). Condition (15a) reduces 
to the inequality 


!/* ; _ I;* / 

2 l 0Cr01 2 UCvll 


5 ^ 01^00 — \hiho 


< 0, (17) 


which can only hold when l m = 0 and |?io| 2 |^n | 2 < 0 . 
Also, (15b) reduces to the inequality 


[L\V}[V,L\ = 


\ho\ 2 0 \ > / c 0 

0 |Z 0 i | 2 / "\0 0 


Thus, to satisfy Condition DS we must select ^oi = 0, 
Zn = 0 and l w ^ 0. As a result, according to Theo¬ 
rem 13, coupling the system with the environment using 
any L of the form of 


» ho i 1 0 , 


will ensure that the state trajectory p t converges to the 
set Zy , which in this example consists of only one ele¬ 


ment p' = 


the operator V. 


, which is the unique ground state of 


We now confirm this funding by directly obtaining the 
equilibrium state of the Schrodinger-picture evolution 
equation (4) under the dissipation control associated 
with the coupling operator L of the form (19) and 
proving its stability. Under L of the form (19), the 
Schrodinger-picture equation (4) becomes 

p = -i [H, p\ + LpL t - ^L*Lp - \pL^L 


\poo \P 01 
~\Pl0 ~\Pll 


\poo —\poi 
\p00 ~\Pll 


— Uio 1 2 Poo ^oozWio — \ poi 9 
hoPooloo ~ \pio9 |Zio| 2 Poo 


^io| 2 Poo 


~ipoi + ^ooPoo^io — \ poi 9 


ipio + ^lofWoo — \pio9 


^io| 2 Poo 


with g = | Zoo 1 2 + | Zoi | 2 - Its equilibrium state p' must 
satisfy the condition p' = 0 , which implies p' 00 = p' 01 = 0 
since Zio ^ 0. As a result, p' 10 = 0 as well. Therefore, 


is the unique equilibrium state under the 


dissipation control of the form (19). It also follows from 
the above calculation that 

Poo.t = e Voo.O) pio.t = Poi,t> 

Pn,t = Pn,o + (1 — e t )poo,o, 

Poi,t = e -( 5 +1)t poi,o + [ e~^ +1)( ' t ^ T) e~ M2r poofldr. 
Jo 

Clearly, p t —> p’ as t —> oo; this implies that p’ is a 
stable equilibrium of the Schrodinger-picture evolution, 
as predicted by Theorem 13. Finally, it is worth noting 
that this equilibrium state is a pure state (since Tr p' 2 = 
1). That is, the designed dissipation control does not 
lead to the decoherence of the system. 

In Example 14, the use of Condition DS also implies 
the satisfaction of (14), and the system is exponentially 
stable. However in general, Condition DS is a weaker 
condition compared to Condition ES, as Condition DS 
does not necessarily lead to exponential convergence. 
This can be illustrated in the following example 

Example 15 Consider a three-level system. Suppose 
we want to engineer the ground-state stability of the 
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Lyapunov operator 


with the generator and the dissipation functional of W 
calculated to be 


V" = 


/ooo\ 
0 1 0 
v 0 0 2 ) 


( 20 ) 


To do this, the coupling operators are chosen as 



^0 10^ 


fo 0 o\ 

L\ — 

0 0 0 

1 L‘2 — 

0 0 1 


yO 0 o) 


l°10 / 


( 21 ) 


Using these values we compute 



o 

o 

o 


(o 0 o\ 

S(V) = 

0 0 0 

, S(U) = 

0 2 0 


T— 1 

1 

o 

o 


0 1/ 


( 22 ) 


Obviously, G(V) < 0. Additionally, we have S(V) > 
\ V. As a result, the Lyapunov operator satisfies the suf¬ 
ficient condition DS. However, (14) does not hold for V 
and any c > 0, hence it cannot be used to establish ex¬ 
ponential stability. In this case, the dissipation will still 
steer the system to the ground state, although the gen¬ 
erator at the first-excited state is zero. In fact, since the 
dissipation strength at the first-excited state is two times 
the dissipation strength at the second-excited state, the 
system will be partly driven to the ground state from the 
first-excited state. However, there is also a possibility 
that the system will be re-excited to the second-excited 
state, which makes the calculation of the convergence 
speed difficult. 

Example 15 provides the evidence that Condition ES 
may not always follow from Condition DS. This example 
shows that asymptotic convergence is not always expo¬ 
nential. 


4 Scalability of the Lyapunov methods 

In this section, we address the question as to how the 
stability of the subsystems scales up when these subsys¬ 
tems interact and are coupled with environments. Asso¬ 
ciated with the subsystems, consider a collection of op¬ 
erators {W\}, W\ y 0, A € A = {1,2,..., N}. Coupling 
between the system and the environments is described 
by coupling operators {L k }, k e A = {1,2,..., A'}. Also, 
consider the operator in the following form 


G(W) = E E L l W ^ - \ L lLkW a - \w x L\L k 

A—1 k 
N K 

= 9 E E(4[WA, L k \ + L{{W X , L k \), (24) 

Z A—1 k 
N N K 

T)(W) = (25) 

A'=lA=l k 

In the space of the subsystems on which W x associates, 
the generator and dissipation functional are 

K 1 1 
G(W X ) = Y L \ W * L * - ~L\L k W x --W x L\L k , (26) 


and 


K 

®(Wx) = Y,i L lw*w*> L k]> 

k =1 


respectively. It follows from (26)-(27) that 


(27) 


N 

G(W) = YS(W x ), (28) 

A=1 

but in general 

N 

S(W)^S(IU A ), (29) 

A—1 


which indicates that the dissipation behaviour may be 
quite different between W and individual W x . 

As noted, the scalability of stability relates to the preser¬ 
vation of stability under the operation of aggregation 
of stable subsystems. In Subsection 4.1 we consider the 
case where each subsystem possesses a Lyapunov ob¬ 
servable that satisfies the conditions derived in the last 
section. Precisely, we consider the situation where the 
coupling operators {L k } are such that each operator W x 
satisfies either Condition ES or Condition DS for this 
subsystem, i.e., each W x is a Lyapunov operator and 
each subsystem, when considered in isolation, asymp¬ 
totically converges to the set of ground states of W x . 
We then derive additional conditions which guarantee 
that W = J2 x =i Wa is a Lyapunov operator for the ag¬ 
gregated system, and hence the entire system is ground- 
state stable. 


N 

w = Y h 0 , 

A=1 


Another way to approach the scalability of the subsys¬ 
tems stability is via studying the total system directly 
using Conditions ES and DS, without imposing the 



ground-state stability requirement on individual sub¬ 
systems and their corresponding operators W \, which 
might be difficult if the system is complex. This case is 
discussed in Subsection 4.2. 

4-1 Scalability of the ground-state stability of each W\ 

Combined with the results from Section 3, the first state¬ 
ment in Lemma 8 can be formulated in terms of the Lya¬ 
punov stability: 

Lemma 16 If Condition ES (respectively, Condition 
DS) holds for each W\ y 0, then the system converges 
to the set of the ground states of W asymptotically. In 
addition, W is a Lyapunov operator. 

Proof. Lemma 16 directly follows from Lemma 8 . 

The above lemma serves as the basis for the scalability 
analysis in this section. The purpose of this analysis is to 
establish conditions for preservation of the ground-state 
stability of W when aggregating ground-stable subsys¬ 
tems. To this end, suppose the ground states of W\ are 
exponentially stabilized by coupling to the fc-th environ¬ 
ment channel; the case of dissipative coupling will be 
considered later. Specifically, suppose the coupling op¬ 
erator L k is such that 

S(W x ) Lk < —c\W\, c x > 0, (30) 

where G{W\)i, k denotes the single-channel component 
of the generator G{X), 

g(X) Lk = L\XL k - l -L\L k X - \xL\L k . (31) 

After the aggregation, the subsystems with which the 
observable W\ is associated, are subjected to other input 
fields and so additional coherent couplings {L k * ,k ^ 
k, k £ A} are induced. Therefore, we need to ensure that 
these additional couplings do not undermine the ground- 
state stability of W\. Formally, this can be achieved by 
ensuring Condition ES is satisfied in the presence of cou¬ 
pling with the environment channels other than k. For 
example, it is sufficient to assume that 

K 

Y S(W x ) Lk , < 0. (32) 

fc'=i,fcVfc 

Clearly, (32) is a sufficient condition to guarantee that 
the satisfaction of Condition ES can be established from 
(30). For this reason, (32) will be referred to as scalability 
condition. This discussion is summarized in the following 
statement. 


Theorem 17 Suppose for each X £ A, there exists an 
L k such that (30) and (32) hold for each Wx- Then, W 
is asymptotically ground-state stable. 

The case where the ground states of Wx are asymptot¬ 
ically stabilized using a dissipative coupling to the k- 
th environment channel can be considered in the same 
manner. In this case, the following statement holds. 

Theorem 18 Suppose for each X £ A, there exists an 
L k such that 

G(Wx) Lk < o, S(WA)L fc > c a IFa (3ca > 0), (33) 

and (32) hold for each Wx- Then, Condition DS is satis¬ 
fied for W and W is asymptotically ground-state stable. 

Proof. In the light of the previous discussion, we 
have GCWx) < 0 by (32). Moreover, T)(Wx)l , = 

[L' k ,,Wx][Wx,L k '} is always non-negative for any 
L k i ,k ^ k, which yields the following relation 

®(WA) > ®{Wx) Lk > cxWx- (34) 

□ 


As Theorem 18 shows, we have dealt with the cross terms 
in T>(W) by introducing a more conservative condition 
(32), which allowed us to engineer the condition (33) 
on each dissipation functional S (II-'a) individually. More 
explicitly, by stabilizing Wx separately and imposing the 
scalability condition (32), we can guarantee the conver¬ 
gence to the set of the ground states of W without using 
the dissipation functional of the total system. 

4-2 Ground-state stability ofW 

As said before, the other approach to the scalability 
problem is to engineer the ground-stability of the total 
system directly. One way to achieve this is by induction, 
by grouping W x , X = 1,... ,n, into W n = YJx=\ w x, 
and considering W n+ 1 as an additional observable. Note 
that we have n £ {1, 2,..., N} = A, and W = Wn- Then 
the algorithm to ensure W is ground-state stable is to it¬ 
eratively achieve for each n the ground-state stability of 
w n+ 1 = nil 1 Wx =W n + W n+ 1 , by synthesizing cou¬ 
pling operators {L k , k = M + 1,..., A'} additional to 
the coupling operators {L k , k = 1,..., M} that ensure 
the ground-state stability of W n . 

Define d n as the smallest eigenvalue of W n , then we have 
W n — d n y 0. Obviously, d-\ = 0. The scalability condi¬ 
tions arising from the above algorithm in the exponential 
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and asymptotic dissipation cases, respectively, are sum¬ 
marized in the following theorems. The proof of these 
Theorems is given in the Appendix. 

Theorem 19 Suppose Q(W n — d n ) < —c(W n — d n ), 
c > 0, is achieved using a set of coupling operators 
{L k , k = 1 W n +i is asymptotically ground- 

state stable if the additional coupling operators { L k , k = 
M + 1,..., A'} satisfy the Lyapunov condition 

K 

g(w n+1 )+ Y S(w n ) Lh 

k=M +1 

< -cW n+1 + c(d n+1 - d n ), c > 0. (35) 


(ii) On the other hand, if 

K 

S(H / A= n +i) + 2 Re([L k , W n ][W\ =n+ 1 ) -^fc]) 

k =1 

> cWa=„+i, c > 0, (40) 

then the conclusion of Theorem 20 holds. 


Proof. Using (38), it follows from condition (39) (re¬ 
spectively, (40)) that (35) (respectively, (37)) holds. The 
statement of the Corollary then follows from Theorem 19 
and 20 . □ 


Theorem 20 Suppose the conditionsQ(W n ) < 0, ®(W n — 
d n ) > c(W n — d n ), c > 0 are achieved using a set of 
coupling operators {Lk, k = 1,... ,M}. The Lyapunov 
conditions to ensure the ground-state stability of W n +i 
are 

K 

Q{W n+ r) + Y S(W n ) Lh < 0, (36) 

k=M+ 1 
K 

2>(W n+1 ) + 2 Y M[L{, W n ][W n+1 ,L k ]) 

k =1 

> cW\= n+ i - c(d n+ 1 - d n ), c > 0. (37) 

Here {Lk, k = M + l,..., K} denote additional coupling 
operators. 

Interestingly, we can further obtain sufficient conditions 
for Theorem 19 and Theorem 20 without knowing the 
value of {d n }. Note that W n+ \ y 0 and d n is the smallest 
eigenvalue of W n . So we have 


If the additional coupling operators satisfy {[L k , W n \ = 
0, k = M + 1,..., K}, then the conditions of Corollary 
21 reduce to conditions (30) and (33). 

In contrast to the scalability approach considered in 
Lemma 16 and Theorem 18, conditions (37) and (40) 
involve the cross-coupling terms 

K 

2YM[LlW n )[W n+1 ,L k }), (41) 

k =1 

and 

K 

2 Y R <[LlW n ][W x =n+h Lk])’ (42) 

fc=l 

These cross-coupling terms show that the condition on 
the dissipation functional of each W\ and the condition 
on the dissipation functional of W do not necessarily im¬ 
ply each other. Therefore, the two methods to achieve 
scalability proposed in this section have different impli¬ 
cations. 


dn +1 — ( 1 ^n+l)p"+ 1 — (Wn + bU„+i) p n+i 
if (fUn)^n-t-l ^ d n , 


We conclude this section with an illustration of the re¬ 
sults of Theorem 19 and the first statement of Corollary 
(38) 21. 


where p™ +1 denotes the ground state of IU„+i- Based on 
(38), we can obtain sufficient conditions which are not 
dependent on {d n }. 

Corollary 21 

(i) Suppose the following condition holds: 


Example 22 Recall the two-level system from Exam¬ 
ple 14, which was concerned with ground-state stabiliza- 

i . 1 o' 

tion of the operator W\ = ^{1 + a Zl ) = 


The Pauli matrices are defined by a x = 



A 0. 


, CTy — 


K 

G(W X 

=n+l) + E g{w n ) Lk < — cWa =n+1 j 

/c—M+l 

(39) 

Then the conclusion of Theorem 19 holds. 


0 -i 
i 0 


1 0 
0 -1 


|. In example 14, the coupling 


operator Li = u_ = 


0 0 
1 0 


was found to satisfy the 


stability condition G(W\)l,i < —W\. Now we consider 
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the extended two-qubit system, on which the operator of 
interests is W = Wi + W 2 with W 2 = |(1 + a Zl a Z2 ) y 0; 
here W\ is the extended operator W\ = W\ <g> I 2 . Also, 
let us extend L\ by letting L\ = L\ ® I 2 . Condition (39) 
with c = 1 in this case takes the form 


G(W 2 ) Ll + £ G(W 2 ) Lk + Y S(Wi) Lh < -W 2 , (43) 


k—2 


k—2 


where L 2 and L 3 are the new coupling operators. Con¬ 
dition (43) can be further simplified 

3 


Y W + < -w 2 - g{W 2 ) Ll 

k—2 

^ 0000 ^ 
0100 
0000 
v o 0 0 lj 


0 0 \ 

0 W' X J 


(44) 


Since W± + W 2 equals 


In this example, since W y 0, the individual W± i2 are 
also asymptotically ground-state stable. In addition, W 
is frustration-free. 


5 Synthesis of the dissipation 

In this section, we introduce the methods to find the 
correct dissipation controls that steer the system to 
the ground states of given candidate Lyapunov opera¬ 
tors. Also, we will show how to calculate the system- 
environment couplings which satisfy the scalability 
conditions derived in Section 4. 

5.1 Synthesis in the case of single dissipation channel 

In this first part, we use a single candidate Lyapunov op¬ 
erator V y 0 and a single system-environment coupling 
operator L as the dissipation control. While we estab¬ 
lished previously that L can be calculated using Condi¬ 
tions ES and DS, these conditions generally lead to non¬ 
linear matrix inequalities. Solving these inequalities for 
large-scale systems is a challenging task. In this section, 
we develop a method to circumvent these difficulties. 


W x +W 2 


1 2 0 0 o\ 
0 10 0 
0 0 0 0 
^0 0 0 ly 


/ 2 0 
\ 0 w' 


letting L 2 , L 3 have the form 


L 2 — 


1 0 

0 a § x3 



W = 


1 ,° ) 


(1 

0 o\ 

0 

0 0 

V° 

01/ 


(45) 


(46) 


reduces inequality (44) to the inequality 

3 

Y0( w ')l* x3 <- w ',w' y 0. (47) 

fc =2 


We introduce a special class of dissipation controls that 
admit factorization L = UV, where U is a unitary oper¬ 
ator. The Reader who is interested in the physical real¬ 
ization of such an operator L can refer to the Appendix. 
In [44] the authors have suggested similar coupling oper¬ 
ators Li \ = UiH\ for the ground-state engineering of a 
Hamiltonian H\. They showed that this class of control 
could form a sufficient condition for ground-state sta¬ 
bility if {Ui} is a set of unitary operators which rotate 
part of the high-energy space with support in H\ into 
the ground-state space, according to [16]. However, it is 
not clear when this rotation exists, and how to solve for 
such unitary rotation. 

In this section we characterize the unitary rotation U re¬ 
quired to establish the ground-state Lyapunov stability. 
Basically, we attempt to solve Conditions ES or DS for 
U. We have 


The easiest way to solve (47) is to further decompose W 
as 


^10 0 ^ 


(0 0 o\ 

0 0 0 

+ 

0 0 0 

^0 0 oj 


(001/ 


(48) 


then the coupling operators are readily computed to be 


L 2 


^0 0 0 0^ 
0 0 0 0 
0 10 0 
v 0 0 0 Oy 


^0 0 0 0 ^ 
0 0 0 0 
0 0 0 1 
v 0 0 0 Oy 


(49) 


g(V) = tfVL - \tfLV - ^VL^L 

= VU^VUV - V 3 (50) 

for single system-environment coupling L. We now con¬ 
sider several special choices for the operator V. 

5.1.1 Special case 1. V is a projection (V 2 = V) 

In many cases, V can be constructed as a projection, 
i.e., V 2 = V. For example, the Hamiltonian discussed 
in [44,30] can be displaced by a constant to generate a 
Lyapunov operator V which is also a projection. More¬ 
over, the two examples of physical relevance considered 


11 



at the end of this section employ Lyapunov operators 
which are projections. 

With the aid of this property, Condition ES can be 
rewritten as 

VUWUV < (1 - c)V, 0 < c < 1, (51) 

which can be regarded as the mathematical formulation 
for the argument in [44]: U should be designed to rotate 
part of the high-energy space into the zero-energy space. 
This shows that our stability results are consistent with 
the physical intuition. 

Now we turn to Condition DS. With L = UV, this con¬ 
dition can be written as 

VUWUV < V, (52) 

and 


where x is an n 2 x 1 vector of free parameters. (V T ®V) + 
denotes the unique Moore-Penrose pseudoinverse [19] of 
V T (g> V. For convenience, we adopt the notations 

i n 2 Xn 2 -(v T ®V) + (V T ®p) = («r 4 ... 4 ) T , 

(V T ® C)+vec(yr^Q) = (&f b\ ... blf, 

(58) 

where the elements {a,i,i = 1,... ,n} are nxn 2 matrices, 
and {bi,i = l,...,n} are n x 1 vectors. According to 
(57), U can be expressed as 

U = (a\x + bi a 2 x + b 2 a n x + b n ), (59) 

which is annxn matrix. The parameters {a.;} and {&,} 
are already known because V, Q are given, and x is de¬ 
termined from the condition U'U = I. The latter con¬ 
dition can be explicitly written as 


©(!/) = [VU\V][V,UV] = -VUWUV + V > cV, (53) 

(3c > 0). 

Obviously, (53) implies (52). More importantly, (53) and 
(51) are the same conditions. As a result, the sufficient 
conditions (14) and (15a) for the ground-state stability 
of V both reduce to the same expression (51) under the 
assumptions L = UV and V 2 = V. In this case, Condi¬ 
tion DS also leads to exponential convergence of (V) Pt 
to 0 . 


/ {ci\X + b^ \ 


{a\X + 6 2 )^ 


(i aix+bi 


a 2 x+b 2 


a n x+b n ) = I. 


\ (a n x + b n y ) 


(60) 

Equation (60) can be further organized as a set of bilin¬ 
ear equations: 


Letting c = 1 in (51) leads to a special case where 
VUV = 0. In particular, a unitary rotation U satisfy¬ 
ing VUV = 0 always exists when stable states are engi¬ 
neered to be the ground states of V [8,44]. 

In addition to the above special case, (51) can be solved 
by making the substitution P = VUV which leads to 
the condition 


X s a\cax + x^a\bi + b\a,iX + b\bi = 1, i = 0,1,... ,n, 
x^ajajX + x^a\bj + b\a,jX + b\bj = 0, i > j. (61) 


The special case where c = 1 and VUV = 0 corresponds 
to b\ = b 2 = ■ ■ ■ = b n = 0. In this case, (61) can be 
simplified as 


(1 - c)V -P f P >0, 0 < c < 1. 


( 54 ) x^a\aiX = l, i = l,...,n, x^a\a,jX = 0 , i> j. (62) 


Since V >z 0 , it ca n be decomposed as V = Q 'Q. There¬ 
fore, P = \/l — <:() is a solution to (54). The synthesis 
problem is transformed to solving 

VUV = 0 < c < 1 (55) 

for a unitary U. (55) is equivalent to 


Example 23 For the purpose of illustrating the differ¬ 
ence between the cases c = 1 and c < 1 , we again con¬ 
sider a quantum two-level system of Example 14 where 
we considered the problem of engineering the ground- 
state stability of the Lyapunov operator (16). With V de¬ 
fined in (16), the Moore-Penrose pseudoinverse of V T (&V 
is calculated to be 


(' V T (g> V)\ec(U) = vec(\/l — cQ), (56) 

vec (U) is the vectorization of an n x n matrix U by 
stacking the columns of U into a single column vector of 
dimension n 2 x 1. The general solution to (56) is given by 


vec(17) = ( V T (g> P) + vec(-\/l — cQ) 

+ (/„ 2xn2 - (V T <g> V)+(V T ® V))x, 


(' V T ® V)+ = 


( \ 0 0 o\ 
0 0 0 0 
0 0 0 0 


\0 0 0 0J 


(63) 


First we solve (62) in the case where c = 1. {a,} can be 


(57) 
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obtained using (63): 


Accordingly, the unitary rotation and coupling are 


a i 


(® oo ®'),*-( o ® io y 

\oiooy \oooiy 


(64) 


With this, (62) can be written as 


U = 


VI - c x 3 


x 2 x 4 


L = 


VT^c O' 
x 2 0 


x 2 = c. 


(71) 

(72) 



^0 0 0 0^ 


^0 0 0 0^ 

X* 

0 10 0 

X = 1, X' 

0 0 0 0 


0 0 0 0 


0 0 10 


yO 0 0 Oy 


v 0 0 0 i ) 


/o 0 0 o\ 


0 0 0 0 

yoio oy 


Equation (72) gives the general form of the coupling 
operator L which satisfies (51) for c < 1. Obviously, (72) 
does not incorporate the special case (68) for c = 1 since 
\x 2 \ < 1. 

5.1.2 Special case 2. V 2 > V 

In this case, the satisfaction of Conditions ES and DS 
still follows from (51) because 

g{v) = vuWuv -V 3 < vuWuv-v < -cv. ( 73 ) 


Parameterizing x as x = (x\ x 2 x 3 x 4 ) T , we arrive 
at a set of bilinear equations 

\x 2 \ 2 = 1, \x 3 \ 2 + \x 4 \ 2 = 1, xlx 2 = 0. (66) 

Particularly, we have x 4 = 0 by (66). Then the unitary 
rotation is 


Additionally, we have 

U f VU < U f V 2 U, (74) 

VUWUV < VUW 2 UV. (75) 

Now it is easy to see that if U satisfies (54), then it also 
satisfies (51), since 


U = ( a 4 x a 2 x) 



M 2 = 1, (67) 


and the desired system-environment coupling L = UV 
is 


L = 



( 68 ) 


Next we consider the case when c < 1. The decomposi¬ 
tion V = Q^Q is not unique, ho wever d ue to the particu¬ 
lar form of (63), (V T ®P) + vec(Vl — cQ) is nonzero only 
if the first entry of Q is nonzero. For example, we can 

choose Q as the square root of V: Q = VV = 

which gives 



Equation (61) transforms to 
13^21 2 = c, \ x 3\ 2 + \x 4 1 2 = 1, 

X* 4 x 2 + y/\ — CX 3 =0, TO > 0. (70) 



VU^VUV < VUW 2 UV < (1 - c)V, c > 0. (76) 

As a result, the unitary solution U obtained from (54) 
could work for both cases V 2 = V and V 2 > V. 

5.2 Synthesis of multiple dissipation channels 

In this section we extend the coupling synthesis approach 
considered in the previous section to construct multiple 
dissipation channels aimed at ground-state stabilization 
of a Lyapunov observable V. We still assume V is a 
projection. Letting L k = U k V , we can re-express (14) as 

K K 

Y, VUlVU k V = V(Y UtVU k )V < (K - c)V, c > 0. 

k =1 k =1 

(77) 

It is easy to verify that Condition DS is still equiva¬ 
lent to Condition ES if we assume the decomposition 
L k = U k V for each L k . This observation leads to the fol¬ 
lowing “no-go theorem” concerning the validity of such 
decomposition. 

Corollary 24 Suppose the Lyapunov observable V is a 
projection. If it satisfies Condition DS but Condition ES 
does not hold, then at least one of the coupling operators 
Lk does not admit decomposition of the form L k = UkV 
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with a unitary Uk ■ In the single channel case, the coupling 
operator L cannot be represented as form L = UV with 
a unitary operator U. 

To illustrate the above result, consider Example 15 
where Condition DS is not equivalent to Condition ES. 
We conclude that the matrix L 2 in Example 15 cannot 
be written as L 2 = U 2 V where U 2 is a unitary operator. 
If this decomposition was possible, we would have 


/o 0 0 N 


fo 0 o\ 

1 Uoo Uq\ \ 

0 0 1 

~ V U w Un ) 

0 1 0 

V° 1 

v 0 0 2 ) 


(78) 


where C/oo is a scalar. It is easy to see from (78) that 

/o 

Un = 2 , Z7oi = [0 0], which implies that U 2 

\10 J 

cannot be a unitary operator. 


The set of unitary operators {Uk} satisfying (77) can 
be calculated if a decomposition such as ( K — c)V = 
Y^k =l QkQk is available. Then {Uk} are obtained by 
solving VUkV = Qk, as did in (55). 


5.3 Scalable dissipations 

In this section we are concerned with a particular class 
of coupling operators of the form Lk = UkW\, where 
W\ y 0, A = 1,2, ...,N and Uk is a unitary operator. 
In this case, there is one-to-one correspondence between 
each Lk and W\ and so we have k = 1,2For 
this type of coupling operators, we can re-express the 
condition [ L k /, W\] = 0, k ^ k, which is one particu¬ 
lar sufficient condition to guarantee satisfaction of the 
scalability condition (32), in terms of unitary operators 
Uk- This leads to a sufficient condition for ground-state 
stability of the operator W = W\, which follows 

from Theorems 17 and 18: 


Corollary 25 Assume \W\,W X '} = 0 for A / A and 
[U k >, W x ] = 0 .If either (30) or (33) holds for each 
W\ 0, then W = Y2\=i Wa is asymptotically ground- 
state stable. 


Proof. The conclusion follows from 

[L k ', TT A ] = [U k >, W X ]W X . = 0, A' ± A. (79) 

Thus G{W\)l , = 0 and so the scalable condition (32) 
used in Theorems 17 and 18 is satisfied. □ 



Fig. 1. The schematic representation of lattices on a torus 
with periodic boundary. The qubits are placed on the edges. 
The qubits 1, 2, 3,4 are connected to a common vertex, and 
the qubits 3,4, 5, 6 form a plaquette. 


If W\ is a projection, the general sufficient condition (32) 
is then expressed as 


W\U^,W\U k 'W\ < W\. (80) 


We now present two examples of application of the scal¬ 
able condition. 

Example 26 Consider the generation of one-dimensional 
cluster state for one-way quantum computation [33]. 
The system is composed of a chain of N qubits with 
nearest neighbor interaction. The one-dimensional clus¬ 
ter state is the ground state of the candidate Lyapunov 
operator defined as 


N-l 

W = E W *’ W * = + !)• ( 81 ) 

A=2 

{ W\ } are commuting due to 

[W x ,W x+ i] = ^ °'zx-i[ (J xxer zx+ i, a zx a xx+1 ]<j zx+2 — 0 . 

(82) 

U\ = ct Za is a solution to W\U\W\ = 0 and (30). Fur¬ 
thermore, we have 


P X ',W X ] = ^[<rz x , ,<rz x - 1 <rxx<rz x+1 ] = 0 (83) 

for A zfi A . By Corollary 25, IT = J2\=2 is 
asymptotically ground-state stable. In particular, 
W can be stabilized to its ground states by select¬ 
ing the following operators as coupling operators 
{L\ = crzx( a zx-i <J nx cr zx+i + !)}• 

Example 27 The toric code [12,25,10] is defined on 
spin lattices. The qubits are placed on the edges, as 


14 


shown in Figure 1. The toric code states can then be de¬ 
fined by the degenerate ground states of the stabilizer 
operators as 

4 6 

a=ik. s= n°v ( 84 ) 

i=l j= 3 

The code states are the ground states of both A and B. 
{ u x } are four X-axis Pauli operators acting on the four 
qubits that connect to one vertex, and {a z .} are four 
Pauli operators acting on the four qubits that form one 
plaquette. Errors can be easily detected and corrected 
using these code states as the computation basis. Nor¬ 
mally, the code states are defined using a large number 
of stabilizer operators, and so the qubits could cover a 
large area of the torus. 

First, we consider two stabilizer operators as out¬ 
lined in (84) and define the Lyapunov operator V 
as V = Vi + V 2 ,Vi = -\a Xl a X 2 a X 3 cr Xi + \,V 2 = 
— ^a Z;t az i cr Z5 a Zg + Performing the similar analy¬ 
sis as in the last example, any of the four operators 
Ui = <J Zi ,i £ {1,2,3,4} can be shown to stabilize V\. 
Particularly, cr Zi ,i £ {1,2,3,4} commutes with any of 
the four operators {cr Zi ,i = 3,4, 5, 6 }. In other words, 
a Zi ,i £ {1,2,3,4} stabilizes V\ without interfering with 
V 2 . By Corollary 25, V) and V 2 are scalable and V is 
asymptotically ground-state stable. 

Suppose the code states are defined by the com¬ 
mon ground states of three stabilizer operators, 
namely, A , B and a third stabilizer operator as 
V 3 = — ^o’ Xl a X7 <J Xs o'x 9 + Vi and V 3 have one common 
edge. In Figure 1, the qubits 7, 8 ,9 should be placed on 
the three edges connecting to the vertex on the left of 
qubit 1. If we use U = a Zl to stabilize V\ , then we have 
Wzi,Vi] = [<? Zl ,-\ <J Xl a X 7 cr Xs a Xg \ ^ 0. The scalable 
condition in Corollary 25 does not hold and so a Zl acts 
non-trivially on V 3 . However, it can be easily seen that 
<j Zl indeed stabilizes V 3 as well, which can be considered 
as a special case where (32) holds. 

6 Conclusion 

We have developed the Lyapunov theory of the ground- 
state stability of quantum systems using a Heisenberg- 
picture approach. This theory is designed to serve as a 
foundation for a theory of stabilization by dissipation, 
which has significant applications in future quantum 
technologies. This theory allows us to engineer the sys¬ 
tems by considering Lyapunov operators and manipulat¬ 
ing inequalities involving the system generators applied 
to these operators, which is a common practice when 
engineering classical and quantum control systems. Sev¬ 
eral issues should be taken into consideration concerning 
practical implementations of this theory. For example, 
the realization of the resulting coupling operators {L^} 


with the available experimental resources is an impor¬ 
tant and interesting problem. Additional Hamiltonian 
control could be introduced if the choices of dissipation 
channels are limited, as was done in [42,43]. The con¬ 
vergence speed is also critical for a large-scale quantum 
system. We have shown that the dissipative stability 
condition DS does not necessarily guarantee exponential 
convergence. Therefore, the scaling of the convergence 
speed as we build up the systems using weak Lyapunov 
conditions will need further investigation. 
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7 Appendix 

7. 1 Proof of Proposition 5 

If X is asymptotically ground-state stable, then (X) Pt —► 
d by Definition 2. Since pt is tight, the limit point al¬ 
ways exists. For any limit point p which is the limit 
of a converging subsequence pt k satisfying p tk —tp, 
we have {X) p > = lim fc ^ 00 (A) ptfe = lim t ^ 00 (A) Pt = d 

and so p is contained in Z x . If pt is converging to Z x , 
then for any converging subsequence pt k of pt we have 
(X) Pt —» (X) p i = d , where p is the limit point. There¬ 
fore, lim t _ >00 (A') Pt exists and equals d. 

7. 2 Proof of Lemma 8 

To prove (i) we observe that if (W\) Pt —> 0, then (W) Pt = 
EJW x ) Pt —>• 0. It remains to show that d = 0. Since 
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W is a finite-dimensional operator, it has finite number 
of eigenvalues. Suppose the smallest eigenvalue of W 
is positive, i.e., d > 0. Then W — d > 0 and ((W — 
d)) Pt > 0 for any state p t . Thus, {W) Pt > d > 0 for any 
state pt (since Trp t = 1). This contradicts {W) Pt —» 0. 
Conversely, if (W) Pf —> 0, then {W\) Pt —* 0 as ( W\) Pt >: 
0 for each A, which proves the ground-state stability of 
each W\. 

7 .3 Proof of Lemma 10 

Since G(V) < 0, (V) Pt < {V) Po and lim t ^. 00 (y) Pt exists. 
Recall that since V is a Lyapunov operator, then Xy — 
{p : (V) p = 0}. Hence Vp 0 € Z v , (V) Pt = (V) Po = 0; 
this implies that Zy is an invariant set. We only need 
to prove that pt will exit the domain {p : {V) p > e} 
for arbitrary e > 0. Now suppose the trajectory p t is 
restricted to a domain {p : (V) p > e} for some e > 0. 
There exists an invariant state pi which is the limit point 
of the tight sequence \ p t > dt [27]. \ p t > dt is the 
mean of pt , and so pj is in {p : (V) p > e}. By assumption, 
(' Q(V)) PI < 0. Let the initial state be the invariant state 
pi. Integrating Q(V) over [0,f] yields 

(V) PI ~ (V) PI = 0 = f{Q(V)) Pl dt. (85) 
Jo 

This leads to a contradiction as (G(V)) Pl < 0. We thus 
conclude that for any e > 0, there exists t(e) such that 

(V) Pt < {V) Pt(e) < e Vi > f(e). That is, V is asymptoti¬ 
cally ground-state stable. 

7 .4 Proof of Lemma 11 

Similar to the proof of Lemma 10 and the proof of Propo¬ 
sition 7 in [27], we only need to prove that for arbitrary 
e > 0, the domain {p : (V) p > e} does not contain in¬ 
variant states pi- Suppose this is not true and there is 
an invariant state pi in the domain {p : (V) p > e}. Con¬ 
sider the positive operator W = V 2 . The generator for 
W is 

g(w) = vg(v) + g(v)v + ®(v). ( 86 ) 

Let the initial state be the invariant state p/. Integrating 
GiW) leads to 

(W) PI - (W) PI = f{vg(v) + g(v)v + ®(v)) PI dt 

Jo 

= 0. (87) 

To establish a contradiction, we use the following iden¬ 
tity 

-vg(v) + -g(v)v + v 2 + \g(v) 2 = [v + -g{v)] 2 . 

c c c z c 

( 88 ) 


Note that by assumption, (V + \g(V)) Pl = (V) PI > 
e, and hence ( [V + yf/(P)] 2 )p / > e 2 > 0 due to the 
positivity of the variance (X 2 ) p — (X ) 2 > 0 for any 
Hermitian operator X. Using ([U + \g(V)] 2 ) Pl > 0 we 
have 

{vg{v) + g{v)v + v{v)) PI 

> {—cV 2 — ^g{V) 2 + D(V)) PI . (89) 

Next, choose a positive number d > 0 such that g(V) + 
d > 0, and then we have —g(V)(g(V) + d) > 0 since 
—g(V) and g(V) + d commute. The latter inequality 
can be written as — G(V) 2 > dg(V). This results in the 
following inequality 

(vg(v) + g(v)v + v(v)) pi 
> (-cv 2 ~- c g(v) 2 + z(v))p I 

> {—cV 2 + ~S{V) + D(V)) PI 
= (- C V 2 + V(V)) PI 

>0. (90) 

The last line of (90) is obtained using the assumption 
cV 2 < Q(V). As a consequence, (87) is not consistent 
with (90). This contradiction shows that for arbitrary 
e > 0, the domain {p : (V) p > e} does not contain 
invariant states pi , hence any trajectory p t must exit 
the set {p : (V) p > e}. This conclusion results in the 
asymptotic ground-state stability of V, which can be 
established using the same argument as in the proof of 
Lemma 10. 

7 .5 Proof of Lemma 12 

Choose a positive number d > 0 such that V — d < 0, 
from which we can conclude V(V — d) < 0 since V and 
(V — d) commute. This can be rewritten as dV > V 2 . 
Thus we have 

V 2 <dV < -D(U). (91) 

c 

By Lemma 11, p* converges to the ground states of V. 

7.6 Proof of Theorem 19 and 20 

For simplicity, first we consider the integration of two 
operators represented by W\ and Wi with IT’i i2 >r 0. 
W = ITj + Wi could be positive definite. Since our Lya¬ 
punov stability results are derived under the assumption 
that the candidate Lyapunove operator has zero eigen¬ 
value, we circumvent this issue by considering the dis¬ 
placed operator W — d y 0, where d > 0 is the small¬ 
est eigenvalue of W. Suppose the Lyapunov condition 
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Q{W\)l 1 < — cW\ has been established using the coher¬ 
ent coupling L\. We are concerned with engineering an 
additional coupling between the environment and the 
part of the system characterized by the observable W 2 
to achieve the following Lyapunov condition for the to¬ 
tal system 

G(W -d) = G(W 1 + W 2 ) < —c(W - d ) 

= -c(W 1 + W 2 )+cd, c> 0. (92) 

Formally, this problem reduces to the that of the synthe¬ 
sis of an additional coupling operator L 2 which couples 
an environment to W 2 . Decomposing (92) yields 

G(W 1 ) Li +G(W 1 )l 2 + J2 S(W 2 ) Lk 

k= 1,2 

< -c(W 1 + W 2 ) + cd, c> 0. (93) 

Since G(Wi)lx < —cWi, a sufficient condition for (93) 
to hold is 

GiW^L, + G(W 2 ) < -cW 2 + cd. (94) 

It follows from the above discussion that W is asymptot¬ 
ically ground-state stable if (94) holds. Similar results 
can be obtained based on the assumption 

G(W!) Ll < 0, S)(Wi) Ll > cW\, c > 0. (95) 

In order to engineer the stability of the combined system 
achieving 

G(W-d)=G(W 1 + W 2 )<0, 

®(W-d) = D(W 1 + W 2 ) 

> c(W -d) = c(Wi + W 2 ) -cd, c> 0, (96) 

we exploit the following relation 

S(E^)= E [LlWx][Wy,L k ] 

AeA fee A;A,A'eA 

= 53®(W a ) + X! [LlW x ][W x ',L k }, (97) 

AeA fee A;A#A' 

to write the second inequality in (96) explicitly as 

S>(W 1 ) + D(W 2 ) + E [L{,W x ][Wy,L k ] 

k, A,A'=1,2;A^A' 

>®(W 1 ) Ll +V(W 2 )+ 53 [L\,W x }[W y ,L k ] 

k, A,A'=1,2;A^A' 

> c(W'i + W 2 ) — cd. (98) 


Accordingly, a sufficient condition to guarantee D(IF — 
d) > c(W - d) is 

®(W 2 ) + 53 [L{,w x \ [W x > , L k \ > cX 2 - cd. 

k, A,A'=1,2;A^A' 

. (99) 

The above methods can be readily extended to consider a 
system that involves an arbitrary number of subsystems 
and observables W\ , resulting in Theorem 19 and 20. 

7.7 Physical implementation of the coupling operator L 

The key component for the proposed method is the abil¬ 
ity to physically realize the engineered coupling L. In 
principle, arbitrary dynamical open quantum system can 
be implemented with high degree of precision given the 
appropriate parameter scaling in the following sense 

lim sup WU^(t)ip — U(t)il>\\ = 0, € H, 

fc->oo 0 <t<T 

(100) 

in which 77 is a subspace of the total Hilbert space. The 
U (t) and U ( t ) are the pre-limit and limit unitary op¬ 
erators given by 

dU {k \t.) = {-iH^dt + dB\t)L (fe) - L^dB(t) 

_ ]-l^L {k) dt)U {k \t), f7 (fc) (0) = I, 

dU(t) = (-iHdt + dB\t)L - DdB(t) 

- Ldt)U{t ),. 17(0) = I. (101) 

The implementation of L is done by adding an ancillary 
qubit to the principal system. The principal system and 
the ancillary qubit are subjected to the following Hamil¬ 
tonian 

H W = kn(La+ + + H s , (102) 

where a + , <j_ are operators of the ancillary qubit. In ad¬ 
dition, the qubit is coupled to an environment via the 
coupling L where k^y is the decay rate 
of the ancillary qubit. In the limit of fast decay of the 
qubit (jfe -4 00 ) the principal system defined on 77 can be 
approximated using L = —~^L, H = Hs- This is 
commonly known in quantum optics as adiabatic elim¬ 
ination, where a fast degree of freedom is eliminated. 
Thus, we have a general approach to synthesize the de¬ 
sired coupling operator L. This method has been sug¬ 
gested for dissipation engineering in [44,26]. A recent ex¬ 
perimental implementation of this method can be found 
in [13]. 
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